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ABSTRACT 

A  model  for  calculating  the  index  of  refraction  of 
atomic  species  present  in  the  cavity  of  a  Xenon  Fluoride 
laser  is  applied  to  Neon  and  Neon* (3s).   The  model  con- 
siders the  variation  from  unity  to  be  a  function  of  the 
absorption  cross  section.   Below  the  ionization  threshold, 
the  cross  section  is  a  set  of  discrete  transitions  between 
the  various  energy  levels  of  the  specie  of  concern.   Above 
the  5p  level,  the  discrete  transitions  become  indistinct 
and  the  continuum  cross  section  function  is  extended  into 
this  region.   The  continuum  cross  section  function  is  a 
polynomial  function  fitted  to  independent  research  results 
The  results  for  Neon  I  agree  with  previous  theoretical  and 
experimental  results.   The  model  is  extended  to  Neon* (3s). 
Comparison  with  the  results  of  independent  research  of  the 
polarizability  and  cavity  phase  shift  derived  from  the 
calculation  indicate  accuracy  is  within  ten  percent- 
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I.   INTRODUCTION 

The  earliest  successful  high  powered  lasers  operated 
in  the  vibration/rotation  transitions  of  C0„  and  produced 
radiation  in  the  infrared,  principally  in  the  10.6  micron 
line.   Extensive  research  and  development  of  technology  in 
the  infrared  laser  field  has  produced  impressive  results 
that  clearly  demonstrate  the  feasibility  of  directed  energy 
weapons  utilizing  these  transitions  for  laser  power  produc- 
tion.  The  relatively  long  wavelength  results  in  large 
optical  components  to  form  and  direct  the  beam.   Even  though 
the  highest  power  lasers  are  of  the  C0„  type,  their  ex- 
tremely large  size  limits  their  usefulness  to  laboratory 
experiments. 

The  shorter  wavelength  of  chemical  DF  lasers  provided 
improved  efficiency  and  reduction  in  size  of  the  power 
production  and  cavity  components.   This  was  the  first  laser 
system  to  demonstrate  target  engagement  and  destruction. 
Still  governed  by  the  laws  of  diffraction,  the  optics  were 
large  and  cumbersome. 

Early  work  on  the  excimer  lasers  had  demonstrated  their 
ability  to  lase  in  the  vacuum  ultraviolet;  however,  problems 
with  mirrors  and  competing  absorption  processes  lead  to  low 
efficiencies  for  these  lasers.   The  channeling  of  atomic 


1.   "International  Defense  Digest,"  International  Defense 
Review,  12:316,  Vol.  3,  1979...  .- 


ionization  and  excitation  into  these  excited  states  was 
found  to  be  very  efficient.   Further  research  by  Velasco  and 
Setzer  produced  interesting  fluorescence  in  rare-gas  halide 
combinations  that  led  Ewing  and  Brau  to  the  demonstration  of 

the  rare-gas  halide  exciplex  laser,  which  is  now  included  in 

2 

the  broad  "excimer"  description. 

The  rare-gas  monohalides  are  simple  diatomic  molecules 
which  are  bound  in  the  excited  state  but  are  unstable  or 
slightly  bound  in  the  ground  state.   The  excited  state  of 
the  excimer  species  can  radiate  in  a  broad  band  that  is  red 
shifted  from  the  wavelength  of  the  parent  excitation.   These 
lasers  have  demonstrated  electrical  efficiency  as  high  as 
ten  percent  in  various  wavelengths,  as  depicted  in  Table  I. 


TABLE  1. 

Rare- 

-gas 

Hal 

ide  Laser 

Em. 

Lssion 

Wavel 

engths 

Peak 

wavelen 

gths  of 

most 

intense 

ban 

d  (nm) 

Ne 

Ar 

Kr 

Xe 

107a 

193 

248 

351,353 

b 

170 

222 

308 

b 

166a,b 

206a 

282 

b 

b 

185a 

b 

252a 

F 

CI 

Br 

I 

(  Has  not  demonstrated  laser  oscillation.) 

(  Predissociates,  hence  emission  is  weak  or  unobservable. ) 


2.   Ewing,  J.  J.,  Rare-Gas  Halide  Lasers,  Physics  Today, 
May  1978,  pp.  32-39. 
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Attempts  to  realize  the  inherent  advantages  of  high 
electrical  efficiency  at  short  wavelength  have  resulted  in 
extensive  research  and  development  activity  in  this  area. 

Extensive  literature  on  the  kinematics  of  the  excitation 

3  4  5 
process  and  associated  flow  effects  is  available.  '  ' 

There  is  a  demonstrated  understanding  of  the  chemical  kin- 
etic processes  involved  in  energy  transfer  and  formation 
of  the  excited  lasing  molecules.   Engineering  development 
has  demonstrated  improved  flow  and  acoustics  effects.   The 
latter  were  essential  to  the  development  of  good  beam  quality 
in  a  demonstration  device.    In  this  case,  good  beam  quality 
is  considered  to  mean  uniform  phase  across  the  cavity  aper- 
ture.  Nonuniform  phase  causes  the  peak  axis  intensity  to 
fall  off  according  to 

^i  =  exp  [-(2tt<J>    /X)2] 
o 

where  <j>     is  obtained  by  squaring  and  averaging  the  phase 
r  ins 

shifts  across  the  aperture.   The  spatial  phase  relation  is 

<f>  =  /  nds  . 


3.  Huestis,  D.  L. ,  Small  Scale  Laser  and  Kinetics 
Program,  paper  presented  at  Rocketdyne  International,  Canoga 
Park,  CA,  15  February  1979. 

4.  Hogge,  H.  D.  and  Crow,  S.  C. ,  Flow  and  Acoustics  in 
Pulsed  Excimer  Lasers,  paper  presented  at  AIAA  Conference  on 
Fluid  Dynamics  of  High  Power  Lasers,  II-4,  Cambridge,  MA, 

2  November  1979. 

5.  Hogge,  H.  Dwight  and  Crow,  Stephen,  Flow  Design 
Concepts  for  Pulsed  Visible  Wavelength  Lasers,  Darpa  Report , 
May,  1978. 

6.  Siegler,  R. ,  et  al . ,  XeF  Demonstrator  Program, 
Rocketdyne  International,  15  February  1979. 
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ELECTRON   BEAM 
PUMP 


OPTICS 

Figure  1.   XeF  Laser  Cavity  Geometry 


,   FLOW 
(Xe:Ne:NF3) 


The  initial  efforts  to  limit  phase  distortion  considered  the 
index  of  refraction  in  terms  of  the  Gladstone-Dale  formula- 
tion 

n  =  1  +  Gp 

where  G  is  the  Gladstone-Dale  constant,  hence  the  concern 

7 
for  flow  and  acoustic  effects.   Hogge  and  Crow   state  that 

flow  uniformity  and  reduced  acoustic  effects  have  achieved 

-5 
Ap/p  =  5x10   ,  the  required  flow  uniformity. 

It  has  been  suggested,  however,  that  the  requirement  for 

flow  uniformity  should  be  stated  in  terms  of  An/n  due  to 

nonuniform  pumping  in  the  cavity.    In  a  typical  flow  mixture, 


7.  Hogge,  ibid. 

8.  Fuhs,  A.,  Etchechury,  J.,  and  Cole,  L.,  Progress 
Report  #2,  XeF  Beam  Quality  Considerations,  submitted  to 
Dr.  J.  A.  Mangano,  DARPA,  2  5  February  1979. 
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the  constituents  of  the  medium  (Ne:Xe:NF3)  are  present  in  a 

9 
fixed  ratio   (750:3:1).    Even  in  a  gas  with  fixed  chemical 

composition,  pressure  and  translational  temperature,  and 

hence,  fixed  Ap/p,  An/n  can  vary  according  to  the  ionic  and 

electronic  excitation  of  the  medium  caused  by  pumping. 

Consider  the  Ladenburg  formulation  of  the  Kramers-Kronig 

,.      .      ,  ,  .   10 
dispersion  relation 


n-1    = 


e2N 


Zirm. 


2 


Z 


N£-    • 


X' 


N 


J      k  kj 


2      fkj 


1    - 


N,  g- 
k°j 


V*. 


1-1 


There  are  £   distinct  species  indicating  specific  electronic 
or  ionic  excitation  present  in  the  cavity.   Index  j  sums 
over  lower  levels  and  k  over  upper  levels  of  possible  con- 
figurations.  Identify  the  mole  fraction  as 


X£   N   ' 


1-2 


the  fraction  of  the  total  molecular  population  represented 
by  species  I.      Finally,  define  a  quantity 


2  ,2 


p2  _     XT.X 


2Trm  c 

e   j   k 


X2   -  X?.  ^ 


1  - 


1-3 


In  terms  of  the  quantities  defined  above 


n  -  1  =  N  Xx/T 


I 


1*1 


1-4 


9.   Siegler,  ibid. 

10.   Mitchell,  A.  C.  G. ,  and  M.  W.  Zemansky,  Resonance 
Radiation  and  Excited  Atoms,  p.  238,  Cambridge,  1934. 
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or  the  index  of  refraction  in  the  cavity  is  the  sum  of  the 
contributions  of  all  the  species  present  in  the  cavity. 
Since  N»  =  NX,, 

n  -  1  =  N,^  +  N2T2+  .  .  .  1-5 

The  partial  density  is  given  by  p^  =  N«m*,  where  m»  is 
the  molecular  weight  of  £th  species;   In  terms  of  partial 
density,  equation  (1-5)  becomes 

h  -  1  =  p1T1/m1  +  p2T2/m2  +. . .      i_6 

Introduce  the  mass  fraction  for  the  gas  mixture 

Y.  =  p»./p,  where  p  is  the  density  of  the  gas  mixture. 

The  result  is 

h :  -  1  =  p[71T1/m1  +  Y2T2/m2  +  ...]         1-7 

Equation  (1-7)  shows  that  n  -  1  is  proportional  to  p; 
however,  equation  (1-7)  further  demonstrates  that  n  -  1 
depends  on  mass  fraction  and  index  of  refraction  of  each 
species.   Even  with  uniform  flow,  the  spatial  phase  relation 
across  the  cavity  will  not  be  uniform.   Accurate  modeling 
of  the  phase  effects  requires  knowledge  of  the  populations 
and  the  contribution  to  the  index  of  refraction  of  the 
species  in  the  cavity.   This  paper  will  explore  a  model 
for  the  index  of  refraction  of  the  atomic  species  in  the 
cavity  and  apply  the  model  to  neon. 

In  Section  IT,  the  formulation  of  the  index  of  refrac- 
tion to  be  applied  in  the  current  work  will  be  derived. 
In  Section  III,  the  concept  of  transition  arrays  and  the 
calculations  of  the  transition  array  oscillator  strengths 
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will  be  reviewed.   The  available  data  for  the  cross  section 
of  Neon  I  (2pb)  and  Neon*(3s)  will  be  reviewed  in  Section 
IV.   The  formulation  derived  in  Section  II  will  be  applied 
to  Neon  I  and  the  results  compared  to  previous  theoretical 
and  experimental  work  in  Section  V.   Section  VI  will  discuss 
the  calculations  for  Neon*(3s)  and  these  results  will  be 
compared  to  independent  experimental  results  in  Section  VII. 
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II.   INDEX  OF  REFRACTION  FROM  DISPERSION  THEORY 

The  ratio  of  the  speed  of  light  in  a  vacuum  to  the 
speed  of  light  in  a  medium  is  known  as  the  index  of  re- 
fraction, 

c 
n  =  v  II-l 

where  v  is  the  phase  velocity.   Consider  a  simple  wave 

E  =  E  ei((ij1:'kx)  II-2 

o 

where  the  phase  relation  is  determined  by 

oot-kx  =   constant  II-3 

Then 

v   =   fe  =   £    .  II-4 

dt        k 

Substitution  of  (4)  into  (1)  yields 

kc  __  _ 

n  =  — .  II-5 

However,  the  index  of  refraction  and  propagation  vector 
are  frequency  dependent ,  hence 

n (fll)  =  MiiL)£  .  II-5A 

go 

In  terms  of  the  relative  permittivity  and  relative 
permeability  of  the  medium 

/fTT    =       / 


For  gas  lasers,  as  for  most  non-ferromagnetic  elements, 
K  does  not  vary  significantly  from  unity.     This  fact 
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yields  Maxwell's  Relation 


n  =  /K   =  /K  (co)  II-7 

e     e 

where  &e  is  the  static  dielectric  constant  of  the 
medium  at  frequency  co. 

The  physical  basis  for  the  frequency  dependence  of  n 
and  the  mechanism  whereby  the  phase  velocity  in  the 
medium  varies  slightly  from  c  lies  in  the  interaction  of 
an  incident  electro-magnetic  wave  with  the  array  of  atoms 
constituting  the  dielectric  material.   In  a  gas,  the 
molecules  are  sufficiently  far  apart  that  each  molecule 
may  be  considered  independent.   The  external  field  sep- 
arates positive  and  negative  charges  which  creates  an 
additional  field  component  known  as  the  electric  polar- 
ization (dipole  moment  per  unit  volume) 

P  =  Nex  1 1-8 

where  N  is  the  number  of  molecules  per  unit  volume  and 
x  is  the  average  relative  displacement  of  charge  due 
to  the  external  field.   The  following  relations  from 
Maxwell's  equations  hold 

D  =  eE  =  £   KPE  =  e  E  +  P  H-9 

o  ®     o 

Simply  stated,  the  electric  displacement  is  proportional 

to  the  external  field.   We  can  write 

p 

Ke  "  1  +  ^~__  11-10 

3  =  /jF~  =   vTL  +  5- 
^e         e  E 
o 


11.   Hecht,  E.  and  Za j ac ,  A.,  Optics,  Addison-Wesley 
Publishing  Company,  1974. 
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where  n  is  the  complex  index  of  refraction. 

a2  -  *  +  A  -  *  +  r#  . 

O  O 

For  more  than  one  dipole  (as  in  array  of  N  atoms) 

5  ■  x  +  I    rr  n-13 

N   ° 

or  the  total  change  in  the  index  of  refraction  is  equal 
to  the  sum  of  the  contributions  of  the  dipoles. 

An  analytical  expression  for  the  average  displacement 
of  the  charges  is  easily  developed  in  classical  terms  of 
the  Lorentz  oscillator  model  of  the  atom.   Although  dev- 
elopment of  this  model  is  properly  done  only  with  quantum 
mechanics,  the  classical  treatment  gives  similar  results. 
In  the  Lorentz  model,  the  outer  electrons  are  considered 

bound  to  a  nuclear  core  by  a  restoring  force  proportional 

2 

to  the  displacement  from  equilibrium  (-m  co   x)«   At 

equilibrium,  the  equation  of  motion  for  an  electron  is 


m 


dx2 


e  dt 


7T     +  m  Y  j*r  +  m  to   x  =  0       11-14 
2     e '  dt     e  o  A 


where  oo   is  the  natural  resonance  of  the  system.   If 
o  J 

the  electrons  are  acted  upon  by  a  harmonic  electromagnetic 

iu)t 

wave  of  the  form  E  =  E  e    ,  the  system  resembles  a 

o     '       * 

classical  forced  oscillator  and  Newton's  second  law  yields 

2 

„,      d  x  dx    ,  2  _      fait  TT   ,  c 

me  ^2   +  me   Y   dt    +  rae%   *  =    eEoe  IJ-15 
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Solutions  to  equation  (11-15)  are  easily  found  by  assuming 
a  solution  of  the  form  of  the  forcing  function.   They  yield 

x  =    2    '2  ,  .    Eoe  n-16 

CO     +  CO    +  1YC0 

o  ' 

Substituting  into  (11-13),  and  canceling  the  field  terms, 

yields 

N£e   /meeQ 

n2    =    1    +X  2,2,.  XI"17 

o  CO  +     CO        +     1YC0 

I  o  ' 

and  we  have  the  form  of  the  complex  index  of  refraction. 
Separating  equation  (11-17)  into  its  real  and  imaginary 
components 

n   =   n   +    iK  11-18 

giV6S  2t_    x,„   2      ,..2 


(N//meQ)(M2-   ^) 
n(co)    =   1+221    7-1 2,     ,      2      2 


K(o>)    -  -  £ 


(co^      -   co    )    +   co    y^ 


2 

^    (N^e    /meQ)(co    y^) 


11-19 


> 


2    Z-    ,       2  2.  2       2 

^       (co         +  co     )    +    (0    Y 


Kramers  and  Kronig  first  derived  comparable  equations 
from  quantum  mechanical  considerations.   Their  equations 

are  found  in  various  texts,  and  the  form  published  by 

12 
Landau  is  shown  for  comparison. 


12.   Landau,  L.  D. ,  and  Lifshitz,  E.  M. ,  Electro- 
dynamics of  Continuous  Media,  translated  from  Russian  by 
J.  B.  Sykes  and  J.  J.  Bell,  p.  259,  Addison-Wesley  Pub- 
lishing Company,  1960. 
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„(„)  M  +  lp[     *fifl  *£  11-20 


(a))    =    _  _  p J       -j  2   da>  11-21 


or      -  a) 


where  P  stands  for  the  Cauchy  principle  value.   Using 
the  fact  that  K(co)  is  an  odd  function 

00 

S(.)  =  1  +  |    f        a'K2M  f  11-22 

J  CO  '  -     CO 

0 
(11-18)  can  be  rewritten  in  a  more  convenient  form 


where 


or 


n(oj)  =  n(a>)  +  iK(co)  11-23 


n(co)  =  1  +  5(a))  11-24 


n(E)  =  1  +  5(E)  11-25 


5(E)  .if    Hin^dE  u_26 

0    E'   -  E^ 


where  the  equation  is  now  considered  in  terms  of  the 
energy  of  the  incident  photons,  E  =  hco ,  and  5(E)  is 
the  variation  of  the  index  of  refraction  from  unity. 


and 


20 


13 
Levinger    demonstrates  the  absorption  coefficient  K 

is  related  to  the  total  cross  section  of  the  individual 

atoms   by 

K(E)    =   jj|5  a(E)  11-27 

or 

2E    (E)    =   Nc"ha(E).  II-27a 

Substituting  into  (11-26) 

Ncli  r°°  cr(E)dE' 


(E)    =  5sl/""  2i 


11-28 


0      E'2    -    E2 

and  we  now  have  an  expression  for  the  variation  from 
unity  of  the  real  part  of  the  index  of  refraction  in 

terms  of  the  total  cross  section  of  the  atom  and  the 

14 
energy  of  the  incident  photon.   Cooper    states  that  the 

majority  contribution  to  the  index  of  refraction  comes 

from  photon  energies  significantly  below  100  Rydberg 

(one  Rydberg  is  13.595  ev),  and  that  in  the  region  from 

zero  to  100  Rydberg,  the  scattering  cross  section  is 

negligible  compared  to  the  absorption  cross  section. 

Up  to  the  ionization  threshold,  the  absorption  cross 

section  is  composed  of  a  series  of  discrete  transitions 

defined  by  the  transition  rules  of  quantum  mechanics. 

Beyond  the  ionization  threshold,  the  absorption  cross 

section  is  a  continuum  to  the  limits  of  integration.   For 


13.  Levinger,  J.  S. ,  letter  dated  4  April  1979  to 
James  Etchechury. 

14.  Cooper,  J.  W. ,  "Photoionization  from  Outer  Atomic 
Subshell.  A  Model  Study,"  Physical  Review,  128:681,  15 
October  1962. 
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a  discrete  transition,  the  total  absorption  cross  section  is 
related  to  the  oscillator  strength  for  the  line  by 


/a(E)dE  =  2TT2e2ti 
/,  .  mc 


line 


nm 


11-29 


The  oscillator  density  for  a  line  is  related  to  the 

cross  section  by 

df 

nm  mr> 

a(E) 


nm 


mc 


dE 


2Tr2e2h 


which  now  defines  the  constant  R  = 


rewrite    (11-28) 


mc 


11-30 


2tt    e   h 


,    and  we   can 


where 


5(E)   =   C 


I 


nm 


E      2-E2 

nm 


+   RP 


/- 

th 


q(E' )dE' 
E'2-   E2 


C   = 


Nch 

ttR 


11-31 


E    is  the  discrete  transition  energy 
nm  OJ 

E'   is  the  variable  of  integration 

and-  energy  units  are  Rydbergs. 

15 
This  formulation  was  used-  by  Migneron    to  calculate 

16 
5(E)  for  Helium  and  by  Liggett.   to -calculate  5(E)  for 


15.   Migneron,  R.  and  J.  S.  Levinger,  "Index  of 
Refraction  and  Sum  Rules  for  Helium,"  Physical  Review, 
13:A646-48,  2  August  1964. 

16..  Liggett,  G.  and.  J.  S.  Levinger,  "Calculation  of 
the  Index  of  Refraction  of  Neon  and  Argon,"  Journal  of 
the  Optical  Society  of  America  ,  58:1:109-113,  January 
1968. 
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Neon  and  Argon.   Their  theoretical  results  agree  well  with 
existing  experiments,  indicating  good  accuracy  for  this 
technique. 
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III.   REVIEW  OF  THE  OSCILLATOR  STRENGTH 

In  order  to  evaluate  the  formula  developed  for  the  index 
of  refraction,  the  oscillator  strengths  in  absorption  for  the 
respective  transition  arrays  are  necessary.   Condon  and 
Shortley  define  a  transition  array  as  the  totality  of  lines 

resulting  from  transitions  between  configurations  (e.g., 

2   2   6      2   2   5  17 

Is  2s  2p  to  Is  2s  2p  3s).   Bethe  and  Salpeter   give  the 

oscillator  strength  in  a  central  field  model  for  a  single 

electron  transition  between  configurations  n  and  m,  defined 

by  quantum  numbers  nl   and  n'£'. 

1   ^max  /T,   „  Nn   2 


f    -  T  oTTT  (E  -E  )R  III-l 

nm   3  21+1    K    n   m  nm 


For  V=l-1,    I        =1,    and  for  l'=l+l,    I        =1+1,    E   and 
'   max   '  '   max     '   n 

E   are  the  energy  eigenvalues  of  the  states  n  and  m. 

%        is  the  radial  matrix  element 
nm 


P  (r)r  P„(r)dr  III-2 


nm  j  n  ^  '   ra 

0 

The  single  electron  oscillator  must  be  multiplied  by  the 
number  of  electrons  in  the  outermost  shell  to  determine 

the  absolute  oscillator  strength. 

18 
Gruzdev    has  calculated  the  radial  transition  integrals 

for  neon  using  a  standard  program  which  calculates  the 


17.  Cited  in  Cooper,  J.  W.  ,  "Photoionizat ion  from  Outer 
Atomic  Subshells.   A  Model  Study,"  Physical  Review,  128:681, 
October  1962. 

18.  Gruzdev,  P.  F.,  and  A.  V.  Loginov,  "Neon,  Radiative 
Lifetimes  for  the  Levels  of  2p5ms,  2p5np,  2p5nd  Qm   =  3-6, 

n  =  3-5)  and  2p54f  Configurations,"  Optics  and  Spectroscopy, 

35:1-3,  July,  1973. 


Hartree-Fock  radial  wave  equation  for  the  configurations 
considered.   These  are  reproduced  in  Table  2. 

TABLE  2:  Transition  Integrals  /  P  (r)rP  (r)dr  for  Neon  Atoms 

&     J      n     m 


Transition 

HF 

Transition 

HF 

3s-2p 

0.43 

3d-2p 

0.11 

3s-3p 

4.76 

3d-3p 

6.3 

3s-4p 

0.35 

3d-4p 

9.22 

3s-5p 

0.13 

3d-5p 

0.11 

4s-2p 

0.17 

4d-2p 

0.081 

4s-3p 

4.25 

4d-3p 

1.6 

4s-4p 

10.91 

4d-4p 

11.4 

4s-5p 

1.07 

4d-5p 

19.2 

5s-2p 

0.10 

5d-2p 

0.06 

5s-3p 

0.98 

5d-3p 

0.82 

5s-4p 

10.3 

5d-4p 

3.22  * 

5s-5p 

19.34 

5d-5p 

17.68 

6s-2p 

0.069 

4f-3d 

10.2 

6s-3p 

0.51 

4f-4d 

15.9 

6s-4p 

12.06 

4f-5d 

1.74 

6s-5p 

18.2 

For  equation  III-l,  an  average  level  found  by  summing 
over  the  levels  of  the  configuration  given  in  Moore  will  be 
used  in  place  of  the  energy  eigenvalues  used  by  Cooper,  i.e., 


2 

E  =  -* 


g.E. 


**i 


III-3 


where  i  indicates  the  various  energy  levels  of  the 
configuration  and  g.  reflects  the  degeneracy  or  quantum 
statistical  weight  of  the  level. 
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While  the  energy  eigenvalues  could  be  found  by  evaluating 
the  radial  Schroedinger  equation,  Ligget's  discussion  of  the 
calculation  of  the  index  of  refraction  cites  difficulty  in 
calculating  energy  eigenvalues.   Use  of  published  informa- 
tion for  E.  was  considered  to  be  more  accurate  and,  in 
addition,  expedited  the  calculations. 
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IV.   THE  ABSORPTION  CROSS  SECTION 

Extensive  research  data  are  available  for  Neon 

2   2   6 

(Is  2s  2p  ).   Since  the  demonstration  of  the  Helium-Neon 

laser,  a  vast  accumulation  of  theoretical  and  experimental 
work  has  occurred.   It  is  safe  to  say  that  the  important 
characteristics  of  the  atom  in  its  ground  state  are  well 
known  to  researchers  in  this  area.   Not  so  well  known  are 
the  characteristics  of  the  metastable  (3s)  configuration; 
as  a  consequence,  the  accuracy  of  the  available  data  for 

the  photoionization  cross  section  is  not  known. 

19 
Ederer  and  Tomboulian   measured  the  total  cross 

section  of  ground  state  neon  from  the  ionization  limit, 

1.58  Rydberg.   Their  results  are  tabulated  in  the  reference 

and  plotted  in  Figs.  2  and  3.   In  his  experiment,  Ederer 

measured  the  transmission  of  ultraviolet  light  with  the 

gas  present  and  with  vacuum.   The  cross  section  is  given 

by  ,        i 

0  =  k  ln  (t^)  IV-1 

where  N  is  the  density  of  atoms  and  L  is  the  path 
length.   To  get  monochromatic  radiation,  he  used  a  normal 
incidence  spectrograph.   To  measure  the  intensity  of  the 
current,  both  a  photographic  plate  and  a  Geiger-Muller 
counter  were  used. 


19.   Ederer,  D.  L. ,  and  Tomboulian,  D.  H. ,  "Photoioniza- 
tion Cross  Section  of  Neon  in  80-600A  Region,"  Physical 
Review,  133 :6A: 1528,  16  March  1964. 
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20 
For  energies  above  12  Rydberg,  Piech  and  Levinger 

have  fitted  the  cross  section  into  a  series  of  power  laws 

The  power  laws  are  fit  to  experimental  data  developed  by 

other  researchers.   The  laws  and  their  respective  limits 

are  listed  below: in  Table  3: 


TABLE  3:   Empirical  Fit  for  Photoionization  Cross  Sections 
Using  Power  Laws. 

Energy  (Ry)  Cross  Section  (Mb) 

13.0-63.7  1.12xlO°E   * 

63.7-227.8  3 . 24xl04E~2 ' 72 

227.8-911.3  1.10xl05E"2,89 


Theoretical  and  experimental  knowledge  of  the  photo- 
ionization of  metastable  rare-gas  atoms  is  rather  scarce. 

21 
McCann  and  Flannery   have  employed  the  single  particle 

central  field  approximation  to  examine  the  photoionization 

of  Ne*(3s).   Their  results  are  plotted  in  Fig.  4.   The 

cross  section  of  the  outer  s  electrons  of  the  metastable 

gas  as  a  function  of  the  ejection  energy  has  been  found  to 

yield  considerable  improvement  over  standard  calculations 

based  on  an  interaction  obtained  from  iteration  of  the 

ground  state  configuration.   These  calculations  demonstrate 


20.  Piech,  K.  R. ,  and  Levinger,  J.  S. ,  "Sum  Rules  for 
Neon  Photoeffect, "  Physical  Review,  135:A332-6,  March  1964. 

21.  McCann,  K.  J.,  and  Flannery,  M.  R. ,  "Photoioniza- 
tion of  Metastable  Rare-Gas  Atoms,"  Applied  Physics  Letters, 
31:9:601,  1  November  1977. 
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that  the  cross  section  is  much  smaller  than  that  of  Neon  I 
and  present  an  optical  window  in  the  2000-3000  A  region. 
Although  not  in  the  immediate  region  of  consideration 
(3510-3530  A),  the  window  is  close  enough  to  have  consider- 
able effect.   Further  inner  shell  ionization  occurs  only 
when  relatively  short  wavelengths  (@  600  A)  are  reached. 
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V.   CALCULATIONS  OF  THE  INDEX  OF 
REFRACTION  OF  NEON  (2p6) 


Equation  (11-31)  was  used  to  evaluate  the  index  of 
refraction  in  the  region  of  interest  (3510-3530  A,  or 
.258  Ry),  at  standard  temperature  and  pressure.   If  the 
cross  section  is  measured  in  Megabarns  (1  Mb=10    cm  ) 
and  the  energies  are  considered  in  Rydbergs  (1  Rydberg= 
13.595  eV),  then  the  respective  value  of  the  constants 
in  equation  (11-31)  are: 

C  =  1.00009xl0~4  Ry2 
R  =  .12895  Ry"1  Mb"1 

The  oscillator  strengths  for  the  2p->(3s,  4s,  5s,  3d) 
transitions  were  calculated  using  Gruzdev ' s  transition 
integrals.   For  the  region  of  interest,  the  contribution 
to  the  index  of  refraction  by  inner  shell  ionization  and 
any  double  excitation  is  negligible.   The  calculations 
agree  well  with  the  calculations  of  Cooper,  except  for 
the  2p->-3d  transition.   However,  this  discrepancy  between 
Cooper  and  Gruzdev  does  not  adversely  affect  the  results 
of  the  calculation  due  to  the  relatively  small  influence 
the  2p-K3d  transition  has  on  the  oscillator  sum  term  of 
equation  (11-31).   The  results  of  the  calculations  are 
given  in  Table  4. 

The  discrete  cross  sections  above  the  5s  or  3d  level 
become  merged  and  insignificant  relative  to  the  numerical 
values  of  the  lower  transitions.   This  is  clearly  demon- 
strated in  the  calculations  of  Table  4  and  in  Figure  5. 
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Figure  5.   Energy  Levels  of  Neon  I  and  the  Region 
of  Discrete  Oscillators  Arrays  (Region 
I)  and  Region  of  Extension  of  the 
Continuum  Integral  (Region  II).  (Scale 
in  thousands  of  cm~l.)   (From  Condon 
and  Short  ley,  p.  302.  ) 
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Table  4.   Transition  Array  Oscillator  Strengths 
for  Ground  State  Neon  (2p^) 


Transition 

Et 

chechury 

2p-*-3s 

0.1632 

4s 

0.0279 

5s 

0.0102 

3d 

0.0239 

Cooper 

0.1632 
0.0264 
0.0090 
0.0372 


Process 

Absorption 
Absorption 
Absorption 
Absorption 


To  account  for  transitions  to  5s  and  above,  the  continuum 
cross  section  data  were  extrapolated  into  this  region;  the 
limit  of  integration  was  adjusted  to  1.51  Ry ,  which  is 
determined  by  energy  of  the  5s  level. 

Adjusting  the  limit  of  integration  in  this  manner  is 
considered  valid,  since  the  discrete  oscillator  strength 
distribution  can  be  considered  as  an  extension  of  the  con- 
tinuum oscillator  strength  distribution.   Cooper  discusses 

23 

this  concept  in  greater  detail. 

In  the  continuum,  a  polynomial  was  fitted  to  Ederer's 
experimental  cross  section  data  in  the  region  1.58   to  3.3 
Ry  and  a  separate  polynomial  function  from  3.3  to  13.0  Ry . 
This  resulted  in  a  set  of  polynomial  functions  that,  when 
separated  into  individual  components,  resembled  the  form  of 
the  power  laws  developed  by  Piech  for  the  region  above  13.0  Ry 


23.   Cooper,  J.  W. ,  ibid 
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The  integral  could  be  represented  as  a  set  of  separate  terms: 

3.3       ,13.0        63.7        227.8         911.8 


q(E)dE'    _ 

E'  ^-E^ 
th  1.51      3.3  13.0        63.7        227.8 


fa^r-f   +  I    +f     +f      +f 


where  the  integral  has  been  terminated  at  the  upper  limit 
of  911.8  Ry. 

From  Figs.  1  and  2  and  Table  4,  it  is  apparent  that  the 
major  contribution  to  the  index  of  refraction  is  made  by  the 
region  near  the  ionization  threshold.   In  the  case  of  Neon, 
the  contribution  to  the  index  from  the  region  1.51  Ry  to 
3.3  Ry,  characterized  by  the  peak  in  the  cross  section  curve, 
is  seventy-three  percent  of  the  continuum  term.   The  region 
3.3  Ry  to  13.0  Ry  contributes  twenty-seven  percent  and  the 
contribution  from  the  region  above  this  is  so  insignificant 
as  to  not  be  considered  in  these  calculations.   This  does 
not  result  in  any  significant  error  and  demonstrates  that 
the  majority  of  the  photoionization  to  the  index  of  refrac- 
tion originates  in  the  region  near  the  ionization  threshold. 
The  small  contribution  at  high  energy  is  illustrated  by  the 
exceedingly  small  cross  section  that  can  be  presumed  by 
extrapolation  of  Fig.  3  to  high  energy.   Comparison  of 
results  in  Table  5  indicates  good  agreement  for  the 
Etchechury  model,  the  more  exhaustive  calculations  of 
Liggett,  and  the  definitive  experimental  work  of  Cuthbertson. 
The  results  demonstrate  that  the  simplification  of  the 
problem  results  in  small  error  and  the  accuracy  of  the 
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Etchechury  model  is  considered  sufficient  for  application 
to  the  electronic  species  Ne*. 


Table  5:   Calculation  of 

6(E)  ■ 

=  n-1 

Author 

Source 

(n-l)xlO4 
0.6705 

E 
0.258 

Ry 

X 
3520 

Etchechury 

Theoretical 

o 

A 

Liggett 

Theoretical 

0.6778 

0.25 

Ry 

3630 

A 

Theoretical 

0.6825 

0.30 

Ry 

3020 

A 

Cuthbertson 

Experimental 

0.6788 

0.25 

Ry 

3664 

A 

Experimental 

0.6812 

0.27 

Ry 

3342 

A 
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VI.   CALCULATION  OF  THE  INDEX  OF  REFRACTION  OF 
NEON*  (3s) 


Equation  (11-31)  was  used  to  evaluate  the  index  of 
refraction  of  metastable  Ne*(3s)  at  .258  Ry  (3520  A)  and 
typical  XeF  cavity  conditions.   The  conditions  were  standard 
temperature  and  pressure  and  an  electron  temperature  of 
13500  K  due  to  pumping  and  kinetic  action  in  the  cavity. 
Using  Boltzmann  statistics  with  the  assumed  electron  tem- 
perature of  13500°K,  the  ratio  N(Ne* )/N(NeI )  is  lxl0~6. 
The  respective  values  of  the  constants  in  equation  (11-31) 

are: 

C  =  1.00009  x  10~10  Ry2 

R  =  0.12895  Ry"1  Mb"1 

Drawing  on  the  previous  section,  the  Etchechury  model 

was  used  to  evaluate  the  index  of  refraction  of  Ne*(3s)  in 

the  region  of  interest  at  standard  temperature  and  pressure. 

The  oscillator  strengths  were  calculated  from  Gruzdev ' s 

integrals  for  the  3s-*-(2p,  3p,  4p,  5p)  transitions,  as  listed 

24 
in  Table  6.   McCann  and  Flannery    indicated  that  inner  shell 

ionization  and  double  excitation  is  negligible. 

The  oscillator  strengths  calculated  by  Cooper's  method 

were  verified  by  means  of  the  polarizability  equation  and 


24.   McCann,  K.  J.  and  M.  R.  Flannery,  "Photoionization 
of  Metastable  Rare-Gas  Atoms,"  Applied  Physics  Letters, 
31:9:599,  November  1977. 
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25 

the  calculations  of  Cohen  and  Schnieder.     They  use: 

2r-2  -_,  f 

5  =  e  h  T  JEE n  VI-1 

a    ra   *-    (E  -E  )*  Vi  x 

em    n  nr 

with  the  results  listed  in  Table  7. 


Table  6.   Transition  Array  Oscillator  Strengths  for  Ne*(3s) 

E   (  Rv  ^ 

Transition    f      nmv  y ;  Process 

3s->-2p      -0.0482    1.230  Emission 

3p       1.080     0.146  Absorption 

4p       0.0106    0.259  Absorption 

5p       0.00129   0.0302  Absorption 


By  comparison,  summing  over  the  oscillator  strength  of 

Table  6  and  adjusting  for  one  electron  value  gives: 

-24   ? 
a  =  32.492  ~    32.5  (xlO    cnT). 

This  agrees  well  with  the  calculations  of  Table  7  and  demon- 
strates the  accuracy  attainable  with  the  Etchechury  model. 

Table  7.   Polarizability  of  Ne* 

-94       ? 
State  (Ne*)   (xlO    cm  ) 

3P2  28.2 

3P  29.2 

3P0  29.2 

1P1  33.4 


25.   Cohen,  S. ,  and  B.  Schneider,  "Ground  and  Excited 
States  of  Ne2  and  Ne+2>"  Journal  of  Chemical  Physics, 
61:8,  p.  3238,  15  October  1974. 


39 


The  contribution  to  the  index  of  refraction  of  the  discrete 

f 
transition  is  6 '(E)  =2_, •   The  summation  term  indicates 

E   2-E2 
that  any  term  having  a  transition  energy  near  the  laser 

cavity  resonance  radiation  frequency  (.258  Ry ,  3520  A)  will  have 
a  significant  contribution  to  the  index  of  refraction.   From 
Table  6,  it  is  apparent  that  the  energy  of  the  3s->-4p  transition 
(E   =  .259  Ry)  is  near  the  cavity  resonance  energy,  and  re- 
quires consideration  of  the  individual  lines  that  comprese  the 
transition  array.  This  implies  that  the  sum  over  the  discrete 
transition  terms  in  equation  (11-31)  is  as  shown  in  Table  8. 


Table  8.  Contribution  of  Discrete  Transitions  to  n-1 

Transition  Terms  in  Series 

3s->-2p  Use  array  term  (negative  for  emission) 

3p  Use  array  term 

4p  Use  summation  over  30  lines 

5p  Use  array  term 

....  Above  5p  included  in  continuum  integral 


In  the  transition  array,  there  are  forty  possible  transi- 
tions between  the  four  3s  levels  (ls«,  Is,,  ls4 ,  lS(-,  in  the 
Paschen  notation)  and  the  ten  4p  levels  (3p.,  where  i=l,10). 
Because  of  selection  rules  on  J,  ten  of  the  lines  are  eliminated, 
leaving  thirty  lines  for  consideration.   The  relative  line 

strengths  for  the  transitions  can  be  determined  from  Table  4  of 

26 
Semenov  and  Strugach.     The  relative  line  strengths  are 


26.   Semenov,  R.  I.,  and  B.  A.  Strugach,  "Determination  of 
Coupling  Coefficients  for  the  np^n's  and  np^n'p  Configurations 
from  Experimental  Data,"  Optics  and  Spectroscopy,  24:,  p.  260, 
April  1968. 
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cated  as  the  upper  number  and  the  corresponding  oscillator 
strength  is  the  lower  number  in  Table  7.  The  levels  of  the  4p 
configuration  are  arrayed  in  order  of  decreasing  energy  across 
the  top.   The  3p5  level  is  at  lower  energy  than  the  3p„. 


Table  9.   Relative  line  strengths  (upper  number)  and  oscillator 
strengths  (lower  number)  for  3s-*-4p  transition. 


E 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

*I0 

3p9 

3p8 

39, 

^6 

3p3 

3p5 

3p2 

3P4 

3px       J 

*2 

.056 
.0015 

0 
J 

0 

0.139 
0.0036 

.833 
.0216 

.111 
.0029 

.667 
.0173 

.333 

.0066 

.833 

.0216 

.333 

4 
.0066 

2s3 

.056 
.0044 

0 
J 

0 
J 

.278 
.0216 

0 
J 

0 

J 

.333 
.0259 

.333 
.0259 

0 

J 

0 

3 

J 

^4 

.250 
.0065 

0 

J 

1.375     .625 
.0356  .0162 

.083 
.0021 

.222 
.0057 

0 

.125 
.0032 

.208 
.0054 

.111 

2 
.0029 

.694     2.333        .292      .097       .750         0  0  .208       .625         0 

2s_  1 

•0108     -0361     .0046  -0016     -0116      J  .0032     .0097      J 

Note:  Paschen  notation  is  used  to  identify  levels. 

The  indices  K,J  refer  to  the  transitions  listed  in  Fig.  7 


In  order  to  clarify  the  contribution  of  the  3s->-4p  transi- 
tion to  the  index  of  refraction,  we  consider  the  summation  term 
in  equation  (11-31)  in  terms  of  cavity  wavelength: 

«'<»  - c  I7T7  =  f4lXJrf2  *nm 

E„„,  -E  A  -A 

nm  nm 

Figure  6  is  a  plot  of  function  VI-2,  using  values  for  f  _ 
ox-  o  nm 

from  Table  9.   Line  width  and  damping  are  not  considered  in 
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these  calculations.   Without  these  considerations,  the  contribu- 
tions to  n  approaches  infinity  at  line  resonance. 

Examination  of  Figure  7  indicates  anomalous  dispersion  at 
351.2,  351.6,  and  351.2.   The  2s.-*-3pg  transition  is  the  source 
at  351.6,  the  2s2-*3p,  transition  at  352.1,  and  the  2s5-*-3p,Q 
transition  at  352.1. 

As  before,  the  discrete  cross  sections  above  5p  level  become 
merged  and  indistinct,  and  the  continuum  cross  section  data  is 
extrapolated  into  this  region.   A  polynomial  was  fitted  to  the 
calculations  of  McCann  and  Flannery  to  generate  a  function  for 
the  integrand  of  equation  11-28  of  the  form  of  the  previous 
section.   The  data  and  the  polynomial  approximation  used  in  the 
region  near  the  peak  cross  section  are  indicated  in  Fig.  4. 
Ionization  threshold  for  Ne*(3s)  is  0.365  Ry  and  the  lower  limit 
of  integration  is  extended  to  0.360  Ry,  determined  by  the  energy 
of  the  5p  level.   The  form  of  the  integral  is: 

00 

f     aCE'  )dE'  z    r'395     +    r3-3 

J  E,2-E2    J  J 

th  -360  -395 

The  integration  extends  only  to  3.3  Ry,  due  to  the  fact  the 
results  of  McCann  and  Flannery  terminate  at  that  energy.   The 
question  arises  as  to  the  extent  of  error  introduced  by  neg- 
lecting the  integral  from  3.3  Ry  to  infinity.   Previous  calcu- 
lations demonstrate  the  order  of  magnitude  of  the  error  intro- 
duced by  termination  of  the  integral  nearer  the  ionization 
level.   Since  data  are  not  available  for  the  cross  section 
beyond  3.3  Ry  level,  an  extrapolation  of  the  relationship 
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"tHii   rftOGkflPl   LhLCJLfi.-:.    THE   CONTRIBUTION    OF    33   TO    4P   TRANSITION    IN    NE- 
THE   WAVELENGTH   OF    THE   ::EF   LASER    RADIATION    IS    352      NANOMETER. 
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Sample  Computer  Output  for  the  Calculation  of 
Contribution  to  the  Index  of  Refraction  of  th 
3s-*4p  Transition. 
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THIS   PROGRAM   CALCULATES   THE   CONTRIBUTION    OF   33   TO   4P    TRANSITION    IN   NE* 

THE    WAVELENGTH    OF    THE    KEF    LASER    RADIATION    IS    351      NANOMETER. 
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Figure  8.      Sample   Computer  Output    for   the   Calculation   of    the 
Contribution   to   the    Index   of   Refraction   of    the 
3s->-4p   Transition. 
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found  in  the  calculations  for  unexcited  neon  was  used.   This 
implies  that  the  cross  section  curves  of  Neon  I  and  Ne*(3s) 
have  similar  shapes  characterized  by  a  peak  near  the  ioniza- 
tion threshold  and  a  sloping  off  to  an  asymptotic  minimum  value 
Intuitively,  this  approach  appears  correct  if  we  compare  the 
shape  of  the  curves  near  the  ionization  threshold.   A  more 
rigorous  discussion  was  not  attempted. 

From  Section  V,  the  contribution  to  the  continuum  integral 
of  the  region  near  the  peak  is  .73  of  the  total  contribution. 
Modifying  equation  (11-31)  by  including  the  weighting  factor 


6(E)  =  C 


nm 


E   2-E2 
nm 


+  BR/ 


3.3 


j(E' )dET 
E,2-E2 


VI-3 


.360 


where  B  =  (l/.73)=  1.37.   Calculation  of  the  variation  of  the 
index  of  refraction  from  unity,  using  equation  (VI-3)  yields 

5(E)  =  -7.79  x  10"10. 
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VII.   COMPARISON  OF  RESULTS  OF  CALCULATION 

The  phase  shifts  across  the  cavity  aperture  per  unit  length 
of  cavity  may  be  obtained  using  the  relation 

£*  =  -«|I2  VII-1 

where  L  and  X  are  expressed  in  consistent  units.   A  lag  in 
phase  relative  to  vacuum  is  treated  as  a  negative  phase  angle. 
Evaluating  equation  (VII-1)  at  3520  A  and  introducing  the  value 
of  5(E)  calculated  in  the  previous  section  yields: 

££  =  2.213xl0"3  m"1 

Accuracy  of  the  calculation  is  limited  by  the  available  cross 
section  data  beyond  the  ionization  threshold.   Improvement  in 
the  calculation  would  result  from  extension  of  the  cross  section 
calculations  used  in  Section  VI  to  an  asymptotic  minimum  value, 

as  has  been  done  in  the  case  of  Neon  I. 

27 
Hunter   ,  in  a  separate  approach  to  the  problem,  calculates 

-5    -1 
the  phase  shift  due  to  Ne*  to  be  2x10    cm 

The  model  developed  in  this  paper  demonstrates  good  agree- 
ment with  independent  research  and  is  applicable  to  any  atomic 
species  present  in  the  XeF  cavity.   The  accuracy  is  dependent 
on  (a)  the  availability  of  cross  section  data  beyond  the  ioni- 
zation threshold  of  the  species  under  consideration,  and  (b) 
estimates  of  the  species  population  in  the  cavity.   Current 


27.   Hunter,  R.  0.,  "Draft  on  Transient  Index  of  Refraction 
in  HgCl  and  XeF  Lasers,"  unpublished  report,  Western  Research 
Corporation,  Los  Angeles,  CA,  p.  17,  April  1978. 
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research  into  development  of  a  populations  model  is  being 
conducted  by  LT  Lonnie  Cole  as  a  continuation  of  this  research 
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